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Spin-orbit coupling (SOC) is central to many 
physical phenomena, including fine structures of 
atomic spectra and quantum topological matters. 
Whereas SOC is in general fixed in a physical 
system, atom-laser interaction provides physicists 
a unique means to create and control synthetic 
SOC for ultracold atoms [it]. Though significant 
experimental progresses have been made 
a bottleneck in current studies is the lack of a 
two-dimensional (2D) synthetic SOC, which is 
crucial for realizing high-dimensional topological 
matters. Here, we report the experimental real¬ 
ization of 2D SOC in ultracold Fermi gases us¬ 

ing three lasers, each of which dresses one atomic 
hyperfine spin state. Through spin injection 
radio-frequency (rf) spectroscopy (^, we probe 
the spin-resolved energy dispersions of dressed 
atoms, and observe a highly controllable Dirac 
point created by the 2D SOC. Our work paves 
the way for exploring high-dimensional topolog¬ 
ical matters in ultracold atoms using Raman 
schemes. 


There have been many theoretical proposals in the lit¬ 
erature for creating multi-dimensional synthetic SOC [1- 
El , so that novel macroscopic quantum phenomena and 
quantum topological states could be studied using ultra¬ 
cold atoms |1514231 . Whereas these proposals have not 
been realized in laboratories, physicists have also just 
begun to explore topolo gical p henomena with ultracold 
atoms in optical lattices |24l427j . Here we use the Raman 
scheme, to produce a highly controllable 2D synthetic 
SOC for an ultracold Fermi gas of Such a SOC al¬ 
lows us to create and manipulate a stable Dirac point on 
a 2D plane, which is detected by spin injection rf spec¬ 
troscopy 0. 

We apply three far detuned lasers propagating on the 
X — y plane to couple three ground hyperfine spin states, 
within the ground electronic manifold, |1) = |F = 


9/2,ruF = 3/2), |2) = |F = 9/2,= 1/2) and |3) = 
\F = 1 {2, mp — 1/2), where (F, mp) are the quantum 
numbers for hyperfine spin states as shown in Fig. la, 
to the electronically excited states. Unlike the tripod 


scheme where a single excited state is considered |9l4l ll 
iMi, in used here, the excited states include a fine- 


structure doublet (^1 liae) and {D 2 line) 

with fine structure splitting of ^ 3.4 nm. Each of two 
D-line components additionally has hyperfine structures. 
New theoretical understandings are therefore required. 

The microscopic Hamiltonian of our system can be 
written as 

3 2 

H = J2i^+ei)\i){i\ + J2Ej\j){j\ 

7 ( 1 ) 

i=l j=l 

where p denotes the momentum of atoms, kj (|kj| = 
2p/\i) and uji are the wave vectors and frequencies of 
three lasers, are the Rabi frequencies, i,j are the in¬ 
dices for the three ground hyperfine spin and the ex¬ 
cited states respectively, si and Pj are the ground and 
excited state energies, n is the total number of the 
excited states and Mij is the matrix element of the 
dipole transition. Different from the proposal in refer¬ 
ence iE^EH, each hyperfine ground spin state here 
is dressed by one and only one laser field, regardless of 
the excited states it is coupled to. A gauge transforma¬ 
tion, 1/ ^ g-2(k 1^)^ can be applied to eliminate 

the phase Oi. All results discussed here are therefore in¬ 
sensitive to the phase difference, and the sophisticated 
and challenging phase lockings are no longer necessary. 

Whereas the standard rotating wave approximation 
gets rid of the time dependence of the Hamiltonian, for 
the far detuned lasers, the excited states can be adia- 
batically eliminated, and the Hamiltonian is written as 
Ha =pl/{.‘^rn) + H^y, 

i=l ^ '' i'^i 

Here, 61 is set as zero (energy reference) for simplifi¬ 
cation, 62 {S 3 ) corresponds to the two-photon Raman 
detuning between Raman laser 1 and 2 (1 and 3), and 
ki = hk.'-. All three Ulip = are real, describing 
the Raman coupling strength between hyperfine ground 
states 1/ \i'). These values could be either derived 

from microscopic calculations (Supplementary Materi¬ 
als), or measured directly in our experiments Q. The 
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FIG. 1: Two-dimensional synthetic spin-orbit cou¬ 
pling. a, Schematic of energy levels of for creating 2 D 
SOC. Each of the three Raman lasers dresses one hyperfine 
spin state from |9/2,3/2) (|1)), |9/2,1/2) (|2)) and |7/2,1/2) 
(13)). Atoms are initially prepared in the free reservoir spin 
state 19/2, 5/2). b, The experimental geometry and the laser 
configuration. The Raman lasers 1 and 2 counter-propagate 
along the y axis and the Raman laser 3 propagates along the x 
axis, which are linearly polarized along z , x and y directions 
respectively, c. Energy-momentum dispersions in the absence 
of SOC ( Qi 2 = = 0). The three parabola are dis¬ 

placed from the origin by ki respectively in the {px,Py) plane, 
e. Energy-momentum dispersions of three dressed states after 
turning on SOC ( —Q 12 = 2 Q 13 = 2^23 = E^), which show 
that the dispersions of two dressed states touch at a Dirac 
point. Here, S 2 = S 3 = 0. d, f are plots enlarging the low 
energy parts of c and e. In e-f , the weights of different hy¬ 
perfine spin states in the dressed ones, the eigenstates of the 
Hamiltonian with SOC, are represented by assigning different 
colors to the hyperfine spin states. Red, blue and green rep¬ 
resent |1), 12), and |3) respectively. The depth of each color 
indicates the percentage of the corresponding spin state in a 
dressed state. 


single-photon recoil momentum kr = 2 T\hj\ and recoil 
energy = k!^!2m are taken as natural momentum and 
energy units. Since the dispersion along the z direction 
is not affected by the lasers, we will focus on the 2D 
Hamiltonian Hxy. 

Using the two lemmas in Supplementary Materials, we 
find out that a doubly degenerate point po exists in the 
momentum space, where po satisfies two independent 
equations 


(ki - kz) ■ Po , ^ ^ 

-h — 02 

m 

m 


^ 12^13 

2^23 

^ 12^23 

2Ui3 


^12^23 

2Ui3 

^ 13^23 

2Ui2 


(3) 


Defining the two dressed states with a touching point 
as a pseudo-spin-1/2, and projecting the Hamiltonian in 
equation m to this pseudo-spin-1/2 near the degenerate 
point Po, an effective Hamiltonian at low energies can be 
obtained straightforwardly. 


HsO — i^^xlPx E ^ylPy)(^x “ 1 “ {^x2Px “ 1 “ ^y2Py)^z ( 4 ) 


where Pi=x,y = Pi —po,i- By rotating the momentum and 
the spin (see Supplementary Material), the Hamiltonian 
can be simplified as 

Hso = KpWx + \p'v^'z^ (5) 

where both A^; and \y are finite. Equation (jS]) describes a 
2D SOC, which is equivalent to the Dresselhaus coupling 
if a simple transformation cr'y^cr'^ is applied. 

Such a 2D SOC directly tells one that the doubly degen¬ 
erate point at Po corresponds to a Dirac point with a 
liner dispersion at low energies. In particular, both the 
amplitude and anisotropy \xl\ can be largely tuned 
(Supplementary Materials). As a demonstration, figures 
le (If) and Ic (Id) show the comparison of the energy- 
momentum dispersion with and without SOC. The latter 
shows that energy dispersions of two dressed states touch 
at a Dirac point. 

In this experiment, a homogeneous bias magnetic field 
Bq = 121.4 G along the 2 ; axis (the gravity direction 
shown in Fig. lb) produces a Zeeman shift to isolate 
these three hyperfine spin states from other ones in the 
Raman transitions, as shown in Fig. la. When three Ra¬ 
man lasers are nearly resonant with three ground states 
( ^^i=i, 2,3 ~ 0), the nearest Raman transitions with other 
hyperfine states are |7/2, 3/2) |2) and |7/2, 3/2) |3), 

which have large two-photon Raman detunings about 

X 60 kHz. Thus we can neglect other hyperfine spin 
states and treat this system as a one with three ground 
nearly degenerate ground states. 

As mentioned before, the Raman coupling strength 
Vtii! includes contributions from all excited states of two 
D-line components as shown in Fig. 2a and can be well 
tuned (Supplementary Materials). When the wavelength 
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FIG. 2 : Controlling the Dirac point by changing the 
wavelength of the Raman lasers, a, Red, green and blue 
curves represent the Raman coupling strength Q 12 , and 
r ^23 as functions of the wavelength A of the Raman lasers. 
Depending on A, ^ 12 ^^ 13^23 can be either positive or negative. 
li is the intensity of the zth laser, b-d, The energy bands for 
the wavelength of 763.5 nm (blue-detuning D2 line), 768.85 
nm (between D1 and D2 line) and 772.4 nm (red-detuning 
D1 line) respectively. The two higher dispersions touch at 
the Dirac point in b and d where ^ 12 ^^ 13^23 is positive. The 
two lower dispersions touch at the Dirac point in c where 
^ 12 ^ 213^23 is negative. 

of the Raman lasers is larger than Di line (770.1 nm) 
or smaller than D 2 line (766.7 nm), fli2l2i3ll23 > 0 (see 
the eigenvalues of dressed states in Supplementary Ma¬ 
terials), and the Dirac point emerges in the highest two 
dressed states as shown in Figs. 2b and 2d. In contrast, 
when the wavelength of the Raman lasers is tuned be¬ 
tween the Di line and D 2 line, ^ 12 (^ 13^23 < 0 and the 
Dirac point emerges when the lowest two dressed states 
become degenerate, as shown in Fig. 2c. This is quite dif¬ 
ferent from the work [l^ , where only one excited state is 
considered and far blue-detuned lasers are used to make 
the degenerate dark states the low-lying ones in the man¬ 
ifold of ground electronic states. 

We investigate the energy-momentum dispersions of 
three dressed states and the Dirac point by spin injec¬ 
tion rf spectroscopy, which uses rf field to drive the atoms 
from a free spin-polarized state into an empty 2 D SOC 
system. We start with a degenerate Fermi gas of 
2 X 10^ at the free reservoir spin state |9/2,5/2) in a 
crossed optical dipole trap. A homogeneous bias mag¬ 
netic field is ramped to Bq = 121.4 G. Then three Raman 
lasers are ramped up in 60 ms from zero to its final value. 
Subsequently, a Gaussian shape pulse of the rf field is ap¬ 
plied for 450 fis to drive atoms from the initial reservoir 
state 19/2, 5/2) to the final empty state with 2 D SOC. rf 
field does not transfer momentum to the atoms and spin 


injection occurs when the frequency of the rf matches 
the energy difference between the initial and final states. 
Since the spin state 19/2, 5/2) is coupled via rf to the state 
| 1 ), rf spectroscopy also measures the weight of the | 1 ) 
state, in addition to the energy dispersions with 2 D SOC. 
Following the spin injection process, the Raman lasers, 
the optical trap and the magnetic field are switched off 
abruptly, and atoms freely expand for 12 ms in a mag¬ 
netic field gradient applied along the x axis. Absorption 
image are taken along the z direction. By counting the 
number of atoms in state | 1 ) as a function of the mo¬ 
mentum and the rf frequency from the absorption image, 
we determine the energy band structure and locate the 
Dirac point. 

Figure 3(a2) shows the momentum-resolved spin- 
injection spectrum when the one-photon detuning of the 
Raman lasers is set at the wavelength of 772.4 nm (red- 
detuning DI) and the corresponding theoretical plot is 
shown in Fig. 3(al). The two higher energy dispersions 
touch at a Dirac point as shown in Fig. 3(a2). To further 
visualize the Diract point, we plot the energy as a func¬ 
tion of Py for various as shown in Fig. 3(a3)-(a5). 
When the wavelength of the Raman lasers is tuned to 
768.85 nm between the Di line and D 2 line, two lower 
energy dispersions touch at a Dirac point, as shown in 
Fig. 3(bl)-(b5). These two wavelengths were used to 
investigate the ID SOC I, i. For the wavelength of 
the Raman lasers at 763.5 nm blue-detuning D2 line, we 
have observed that two higher energy dispersions touch 
at a Dirac point, similar to the 772.4 nm case. We also 
perform numerical calculations for the eigenvalues of the 
Hamiltonian in equation ([2]) and (|5j) according to experi¬ 
mental parameters and have found out a good agreement 
between theory and experiments. 

The Dirac point created in this scheme is robust, in 
the sense that it moves in the momentum space with¬ 
out opening a gap when experimental parameters change. 
This comes from the fact that the Hamiltonian in equa¬ 
tion (| 2 ]) is real, and consequently, the SOC in equation (|5]) 
cannot contain the ay term. By modulating Si, Cli and k^, 
an arbitrary trajectory of the Dirac point in the momen¬ 
tum space can be designed in principle. As a demonstra¬ 
tion, we fix the Raman detuning S 2 and measure the posi¬ 
tions of the Dirac point on the {px^Py) plane as a function 
of S 3 , as shown in Fig. 4a (blue dot). When fixing ^3 and 
changing S 2 , the Dirac point move along a different line. 
We also measure the energy at the Dirac point, as well as 
the energy separation AE from highest dressed state, as 
a function of the Raman coupling strength. Zero energy 
in Fig. 4b corresponds to the Zeeman energy splitting be¬ 
tween |9/2,5/2) and |9/2,3/2). The energy of the Dirac 
point decreases and AE increases linearly with increasing 
the Raman coupling strength (inset of Fig. 4b). 

In this work, we have realized a highly controllable 
2 D SOC and a stable Dirac point for ultracold fermionic 
atoms. In the future, we will develop technologies to 
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FIG. 3: The energy dispersions of dressed atoms 

measured by rf spin-injection spectroscopy. al-a7, The 

one-photon detuning of the Raman lasers is set at the wave¬ 
length of 772.4 nm (red-detuning Dl). Q 12 = 3.58F^r, = 

-3.94^;^, Q 23 = -4.66^;^^, 62 = -b.UEr, = -3.23^;^. bl- 
b7, The wavelength is set at 768.85 nm between the Di line 
and D 2 line. Q 12 = —4.96F^r, f^i3 = 5.46F^r, f^23 = 6.46F^r, 
82 — —5.2Er, S3 = —2.13Er. al and bl are the theoretical re¬ 
sults using the realistic experimental parameters. a2 and b2 
are experimental results measured by rf spin-injection spec¬ 
troscopy. a3-a5, The cross-section drawings of (al) and (a2) 
in the energy-p^ coordinates for different quasimomentum px 
Red dots and solid curves are experimental and theoretical 
results respectively. The inset in a4 is plot enlarging at Dirac 
point. Here, A^i = 0, Xx 2 = —0A2kr/m, Xyi = 1.5krlm, and 
Xy 2 = —0.07kr/m are obtained by the experimental param¬ 
eters. b3-b5, The cross-section drawings of (bl) and (b2). 
The inset in b4 is plot enlarging at Dirac point. Xx,i, Xx ,2 
Xy^i and Xy ^2 are the same as those in a, since ^4i2/f^i3 and 
^ 12/^23 remain the same. 



[unit ofkj 



FIG. 4: Moving the Dirac point, a, controlling the Ra¬ 
man detuning S 2 and 83 ^ the Dirac point moves along certain 
trajectories in the momentum space. Blue dots: 82 is fixed as 
—5.45F^r and 83 is tuning at the range of [—2Er^ 0.45F^r]. Red 
dots: 83 is fixed as —0.2^^^ and 82 is tuning at the range of 
[—7F^r, — 2.5F^r]- The other parameters are D 12 = —4.96F^r, 
Qi 3 = 5.46F^r, f ^23 = 6.46F^r. b, the energy at the Dirac point 
and the corresponding energy of the highest dressed state with 
the same momentum po as a function of the Raman coupling 
strength. —D 12 = 0.91^13 = 0.77^423 — 82 = —3.2Er and 

83 = — 2.13F^r- The left inset in b shows the energy disper¬ 
sions, The right inset in b represents the energy separation 
AE between the Dirac point and the highest dressed state as 
a function of the Raman coupling strength. Here, the one- 
photon detuning of three Raman lasers is at the wavelength 
of 768.85 nm. The error bars represent the standard deviation 
of repeated measurements. 


open and control the gap at the Dirac point to gen¬ 
erate an effective perpendicular Zeeman field, so that 
interesting topological and other exotic superfluids will 
rise from s-wave interaction (2 qI-[ 23|. The heating effect 
due to spontaneous emission induced by Raman lasers 
is neglected here, since we inject atoms to states with 
SOC that are initially empty. Such a technical problem 
will become crucial if atoms are initially prepared in the 
states with 2D SOC. It may be overcome by using atomic 
species with the absence of the spontaneous emission, or 
by exploring the quantum quenching techniques so that 
topological properties may be extracted from quantum 
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dynamics in a relatively short time scale [1^ . 

Method: 

Experimental set-up. After sympathetic cooling 
with bosonic in a quadruple-Ioffe configuration 

magnetic trap, in the 19/2, 9/2) state and ^^Rb atoms 
in the |2,2) are transferred into an optical dipole trap 
formed by two orthogonal 1064 nm laser beams. The 
Fermi gas is further evaporatively cooled to TjTp ^ 0.3 
about 2 X 10^ with ^^Rb atoms in the optical trap 
[ 2 ^, where Tp is the Fermi temperature defined by 
Tp = and the geometric mean of trap¬ 
ping frequencies ^ 27r x 80 Hz in our system, N is the 

number of fermions. After removal of ^^Rb atoms, the 
fermionic atoms are transferred into the state 19/2, 5/2) 
via a rapid adiabatic passage induced by a rf field of 80 
ms at 19.6 G. Here, the transition is addressed by a rf 
ramp that starts at 6.56 MHz to the end of 6.28 MHz. 
Then a homogeneous bias magnetic field along the z axis 
(gravity direction) is ramped to Hq = 121.4 G by a pair 
of coils operating in the Helmholtz configuration. 

Raman lasers. Three Raman lasers are derived from 
a GW Ti-sapphire single frequency laser. Two Raman 
beams 1 and 2 are frequency-shifted around +201.144 x 2 
MHz and +220.531 x 2 MHz by two double-pass acousto¬ 
optic modulators (AOM), respectively. The Raman 
laser 3 is sent through two AOMs with double-pass and 
frequency-shifted —212.975 x 4 MHz. Then three Ra¬ 
man beams are coupled into three polarization maintain¬ 
ing single-mode fibers respectively in order to improve 
stability of the beam pointing and achieve better beam- 
profile quality. Behind the fibers, these three Raman 
beams have the maximum power 80 mW for each beam 
and they overlap in the atomic cloud with 1/e^ radii of 
200 fim. The Raman lasers 1 and 2 counter-propagate 
along the y axis and the Raman laser 3 propagates along 
the X axis, which are linearly polarized along z , x and y 
directions respectively, corresponding to tt, cr and cr po¬ 
larization relative to the quantization axis z as shown in 
Fig. lb. 
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SUPPLEMENTARY MATERIAL 


2D synthetic SOC and the Dirac point 


Despite its simple form, equation ( 2 ) in the main text leads to significant results of a highly controllable Dirac point 
and a 2 D synthetic SOC. To see this fact, we first present two lemmas. 


Lemma 7, a real 3x3 matrix 


_ 1 ^ 121^13 _ 1^12 _ 1^13 \ 

2023 2 2 

_ Ol 2 _ O12O23 O23 

2 20i3 2 

_ Ol 3 O23 O13O23 

2 2 2 O 12 


(6) 


has two degenerate eigenstates |A)and \B) with zero eigenenergy Ea = E'b = 0, 

1 ^) = ^( - f^ 23 |l) + fil 3 | 2 )) 

IB) = + 


(7) 


where Ni = + Q 23 ^2 = ^13 + ^23 + normalization factors. The third eigenstate 

\C) has the energy 

1 


Ec = - 


2Qi2Qi3^23 


(Qi2^13 + ^12^23 + 


2 ^2 \ 
13^^231 


(8) 


For positive and negative ^ 12 ^ 13 ^ 23 , \C) is the ground or excited state respectively. 
Lemma 2^ there always exist a finite po to satisfy a set of three equations 


(po - kj)^ 

2m 


Si =e- 


^ii' ^ii 


i = l,2,3 


(9) 


where e is a constant, i ^ i' ^ i", provided that Liu/ ^ 0 for any i,i', and none of these three momenta = fik- is 
parallel to other ones so that the moentas ki — k 2 , k 2 — k 3 , and ki — k 3 , two of which are independent, could span a 
plane. 

These three equations are equivalent to 


Po /I 1 \ I /X X \ ^12^13 , ^12^23 

-ki - k 2 + (5i - 82 ) = —— 

m ^^23 

Po /I 1 \ , /X X \ ^12^23 , ^13^23 

-k2 - k3 + ^2 - ^^3 =-- 

m 2 Ui 2 

Po /I 1 \ I /X X \ ^13^23 , ^12^13 

' (^0 - ki) + {S 3 - Si) = -—-+ —-, 

m zUi 2 ^^23 


( 10 ) 

( 11 ) 


( 12 ) 
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only two of which are independent, since the summation of all three become zero on both sides of the equation. Thus 
it is sufficient to consider only the first two equations. 

If the two momenta ki — k 2 and k 2 — ks are not parallel to each other, they span a plane. Define this plane as the 
x — y plane, the z-component of po drops off from these equations, and we could consider a two-dimensional problem. 
The solutions of equations [10] and im form two lines on this plane. The intersection then uniquely determines the 
location of po on this plane, and the corresponding energy of this degenerate point is e. Note that regardless of the 
microscopic parameters, there is always a trivial solution of po at infinity, which is irrelevant to our discussions here. 

Project the Hamiltonian to the subspace of 1^4) and |H), the effective Hamiltonian is, 

_ ( {A\H\A) / 7^(^^23<^ip + i^i3«2p) 73^(e2p-eip)ili3Si23 \ 

' - I, (B|JJ|2!) (B|ff |B) ) - - 3,p)n,3n23 (ni.^ip + nfjejp + / 

where e^p = (p — ki)^/(2m) Expand near the degenerate point po, and redefine p = p —po, and k^ = k^ — po, 
then the effective Hamiltonian can be written as: 




mN‘-, 




(fJlgki +fjf3k2) -p 
f^l.'if223(kl — k2) • p 


;^ 13 ^ 23 (kl - k 2 ) • P 


^iV| 


(^23^1 + ^I3k2 + 


^12 


13+^23 

^13^23 


O'P. 


e • /2 


x2 


(14) 


where 12 x 2 is the identity matrix. 

Dropping the terms that are independent on spin. He can be in general written as 

He — {^xlPx H“ ^ylPy)^x {^x2Px ^y2Py)(^z 

where the coefficients A^^i, Xx 2 ^ A^i, Xy 2 depend on the microscopic parameters such as k^ and ftu'. 

To simplify the notations, we chose a coordinate so that kiy = k 2 y. Then the effective Hamiltonian ([T4|) can be 
written as 

He = apxCTx + bpxCTz + cpya^ + el (16) 


where a, h, c and e are all constant: 

1 


a = 


mNiN2 


f^l3^23 {klx — k2x) 


^ “ 2 mAf 2 (^ 23 ^ 1 * + 2 mAf| (^ 23 ^ 1 ^ + ^^ 3^202 + ^ 302 ^ 


c = 


2 mlV 2 + ^?3fc2y) 2 mlV 2 (^i3^1J/ + ^?3fc2j; + ^3j;) - (^U/ ^3j;) 


When (ki — k 2 ) x (ki — ks) 7 ^ 0, a 7 ^ 0 and c 7 ^ 0 are satisfied. 

We rotate p about the Pz axis by p and the spin about the (jy axis hy 0. If p and 0 satisfy 

(a^ + 6 ^ — c^) sm 2 p + 2 bccos 2 p = 0 , (a^ — b^— c^) sin46> -f- 2a6cos46> = 0 , 


the effective Hamiltonian is simplified as 


Hso — ^xPx^x ffi ^yP^er 


y^z^ 


(17) 


(18) 


where 


A^ = acospcos20 — bcospsm20csmpsm20 

Xy = asmpsm20 ^ bsmpcos20ccospcos20 (19) 

A special case is a perfect symmetric configuration. Si = 62 = Ss^ fti 2 = fli 3 = ^^23 so that k^ form an equilateral 
triangle. The Hamiltonian is Hso = 2 v^m ^ Py^z) where d = \ki — kj|. 
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Calculation of Raman coupling strength 

The Raman coupling strength can be expressed as [l| 


^ij — 


/q I evq \ F ", rrip) {F ", \ evg \ ) 

JFceo ^ A 


( 20 ) 


where, Iq = yjli • / 2 , and the Ii is the intensity of each Raman laser light, c is the speed of light, eo is the permittivity 
of vacuum, e is the elementary charge, and q is an index labeling the component of r in the spherical basis. A is 
one-photon detuning of Raman lasers. and are two ground hyperfine spin states coupled by a pair 

of Raman laser. is the middle excited hyperfine spin state in the Raman process. 

Here, we consider atoms, the excited states of the laser-atom coupling involve a fine-structure doublet 4^Pi/2 
(Pi line) and {D 2 line) with fine structure splitting of ^ 3.4 nm as shown in Fig. O The Raman coupling 

strength can be written with two contributions from Pi line and P 2 line 0 , 


^ij — 








Ad2 




( 21 ) 


F'',m'i 


Here, Am and Ad 2 are one-photon detuning of Raman lasers relative to the Pi and P 2 line respectively. It is 
required in here that Am and Am are very larger than the hyperfine energy splitting of the excited state. 

The Raman coupling strength for Pi and P 2 line can be simplified as 


Jq (FJ , ml I erg \ F ", m'p ){F'', m'p \ eVg \ , m^p) 

h?ceo A/ji 


Y 

F • S' 


h?ceQ Am 


(J = l/2|er,|J" = l/2)2, 


( 22 ) 


^‘‘ij 


Iq {F^,m^p\ erg \ F ", m'p ){F", m'p \ erg \ F ', m'p) 


h?ceo 


lo Q 




JFceo 


^D2 


Ad2 

-{J = l/2\erq\r = 3/2f. 


(23) 


Here, is the hyperfine dipole matrix element between the ground and the excited hyperfine state 

|P",m^) depending on tt or transition, and (J = l/2|erg|J" = 1/2) and (J = l/2|erg|J" = 3/2) are the PI and 
P2 transition dipole matrix element respectively, whose values all can be found in Q. Therefore, the Raman coupling 
strength ( 2 i 2 , ^13 and r 223 as a function of one-photon detuning can be obtained from summing contributions from 
all the excited hyperfine states of two D-line components using Eqs. [2TII231 as shown in Fig. 2a in the main text. 
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F=5/2 



FIG. 5: (Color online). Schematic of energy levels of 






















